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Abstract. Evidence for excess extreme ultraviolet 
(EUV) emission over a tail of X-ray gas bremsstrahlung 
emission has been building up recently, but in some cases 
remains controversial, mostly due to the moderate quality 
of the EUV data. In order to improve the signal to noise 
ratio in the EUV, we have performed the wavelet anal- 
ysis and image reconstructions for five clusters of galax- 
ies observed both at EUV and X-ray energies with the 
EUVE and ROSAT satellites respectively. The profiles 
fS| ' of the EUV and X-ray reconstructed images all differ at 
. a very large confidence level and an EUV excess over a 
thermal bremsstrahlung tail is detected in all five clusters 
(Abell 1795, Abell 2199, Abell 4059, Coma and Virgo) up 
to large radii. These results, coupled with recent XMM- 
Newton observations, suggest that the EUV excess is prob- 
ably non thermal in origin. 

fSJ ' Key words: galaxies: clusters: individual (Abell 1795, 
Abell 2199, Abell 4059, Virgo, Coma) — Radiation mecha- 
CO ' nisms: general — X-rays: galaxies: clusters — Cosmology: 
. observations 
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ASTRONOMY 

AND 
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1. Introduction 

The Extreme Ultraviolet Explorer (EUVE) has detected 
emissiqn from a low dusters of galaxies In Ilk! ~ 70-200 



clusters could exclude the presence of warm gas at a 2a 
level (Dixon et al. 2001). On the other hand, non ther- 



cV energy range. By order of discovery, the se were: Virgrj 



Lieu e * L al. 1996l | ; | Borghofor oL al. 2000a|) Coma flEe 



et al. 1996a|), Abell 1795 flMittaz et al. 199S|), Abell 2199 



mal models ( Barazin & Lieu 1998; Bowyer fc Berghdfer 



1998; Lieu et al. 1999a; Ensslin et al. 1999; At 



Volk 2000| ; |Brunctti et al. 200 1| ) seem able to account 



toyan 



flBowyer et al. 1998j ), Abell 4059 flBerghofer et al. 2000bD 
and Fornax (Bowyer et al. 2001). The problem of the phys- 
ical origin of the EUV emission was immediately raised. 
Excess EUV emission relative to the extrapolation of the 
X-ray emission to the EUVE energy range was detected in 
several clusters, suggesting that thermal bremsstrahlung 
from the hot (~ 10 8 K) gas could not be entirely responsi- 
ble for the EUV emission. Although there has been some 
controversy on the reality of this EUV emission, it now 
seems well established for the following clusters: Coma 



(Lieu et al. 1996a|), Virgo (Lieu et al. 1996b|; Berghdfer 


etal. 2000a 


Bonamente et al. 2001a 


), Abell 1795 (Mittaz| 


etal. 1998; 


Bonamente et al. 2001a), Abell 2199 (Bowyer 


et al. 1998; 


Kaastra et al. 1999; 


Lieu et al. 1999b) and 


Abell 4038 ( 


Bowyer et al. 1998). 



This excess EUV emission can be interpreted as due to 
two different mechanisms: thermal radiation from a warm 
(10 5 — 10 6 K) gas, as first suggested by Lieu et al. (1996a), 
or inverse Compton emission of relativistic electrons either 
on the cosmic microwave background or on stellar light 
originating in galaxies (Hwang 1997; Bowyer & Berghdfer 
1998; Ensslin & Biermann 1998; Sarazin & Lieu 1998; En- 
sslin et al. 1999), or both. Note that although the existence 



of a multiphase intraclustcr medium is plausible ( Lieu et 



al. 200C ; Bonamente et al. 2001b), there are several dif- 



ficulties with the thermal model. The most serious one 
is that since gas in the temperature range 10 5 — 10 6 K 
cools very rapidly, a source of heating for this gas is neces- 
sary, and only one model has been proposed to account for 
this gas heating (Fabian 1997). Besides, the non-detection 
by the Far Ultraviolet Spectroscopic Explorer of the OVI 
AA1032,1038 lines in the spectra of the Coma and Virgo 
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for the excess EUV emission, but require very high cosmic 
ray pressure. 

Following a controversy that the extended EUV emis- 
sion from some clusters could be simply an artefact of 
bad background subtraction (Bowyer et al. 1999; Lieu et 
al. 1999b) we have performed a wavelet analysis and re- 
construction of the EUVE images of several clusters. Our 
aims were to confirm the detection of EUV emission as 
far as possible from the cluster center, to derive accu- 
rate EUV and X-ray profiles and to confirm the existence 
and radial distributions of the soft excesses over thermal 
bremsstrahlung in these clusters. 

We present here the wavelet reconstruction of the im- 
ages of these five clusters, both in the EUV (EUVE im- 
ages) and in X-rays (ROSAT PSPC images). Abell 4038 
was also analyzed but was discarded due to its low sig- 
nal to noise ratio in the EUV. We then compare the EUV 
and X-ray profiles of the images thus obtained, arguing 
that if both emissions have the same physical origin it 
should not be possible to distinguish the shapes of these 
profiles. In fact, we will see that they are found to be dif- 
ferent in all clusters, and that soft excesses over thermal 
bremsstrahlung are detected in all five clusters with vari- 
ous radial distributions. 

The data and method are presented in Sect. 2. Results 
are presented in Sect. 3 and discussed in Sect. 4. 



2. The data and image analysis 

2.1. The data 

Five clusters were observed in the EUV range with the 
EUVE satellite and in X-rays with the ROSAT PSPC 
satellite (archive data). The exposure times and main clus- 
ter characteristics are given in Table 0. 
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Cluster 


EUVE total 


ROSAT exp. 


ROSAT 


Redshift 


kT x 


Scaling 


Cooling 




exp. time (s) 


time (s) 


Obsld 




(keV) 


(kpc/superpx) 


flow 


A 1795 


158689 


33921 


rp800105 


0.063 


5.9 


32.7 


yes 


A 2199 


93721 


34633 


rp800644 


0.0299 


4.1 


15.8 


yes 


A 4059 


145389 


5225 


rp800175 


0.0478 


4.0 


25.0 


yes 


Coma 


60822 


20112 


rp800005 


0.023 


8.7 


12.2 


no 


Virgo 


146204 


9135 


rp800187 


0.003 


2.4 


1.6 


weak 



Fig. 1. Virgo long exposure - arrangement of boxes in 
order to obtain the image ready for the wavelet recon- 
struction. Note that the different boxes here are only il- 
lustrative of the algorithm we used. 



Fig. 2. Abell 1795 arrangement of boxes in order to obtain 
the image prior the wavelet reconstruction in the case of 
a co-added image. The central circle has a radius of 16 
arcmin, which we will take as the limit of our analysis. 



In order to have comparable spatial resolutions in the 
EUV and X-rays, the EUV images were rebinned 4x4, 
leading to a "superpixel" size of 0.3077 arcminutes (18.5 
arcseconds) . X-ray images were also rebinned to the same 
18.5 arcsec pixel size in order to allow a direct comparison 
of the EUV and X-ray profiles. 

The linear scale per superpixel at the cluster distance, 
estimated with Ho = 50 km s _1 Mpc -1 and qo=0 is given 
in Col 7 of Table |i] 



Fig. 3. Enlarged and rebinned image of Abell 1795. The 
total size of the image is 140x140 arcmin 2 . 



perimposed on the raw image. The largest one, denoted 
by (A,B,C,D) follows the outer detector shape closely. Its 
edges are detached from the detector boundaries to avoid 
the steepening of the background near the edge, which 



would bias the noise estimation (see the figure 2 of Lieu et 



al. 1999c). In this box, we define two rectangles (denoted 



2.2. The EUVE data 

The EUVE/Deep Survey photometer (DS) images were 
extracted from the raw event lists with the customary lin- 
ear bin size of 13 pix/arcmin. The rectangular shape of 
the Lex/B filter (~70-200 eV) results in images where the 
length of one side of the image greatly exceeds the other 
side (see Fig. 1). When more than one observation was 
available for each target, images were coadded to improve 
the S/N; the rectangular shape often results in only partial 
overlap between the exposures leading to diamond-shaped 
geometries (see Fig. 2). Without taking into account any 
edge corrections, this peculiar shape of the EUVE images 
may be a problem for any analysis involving an isotropic 
distribution of the data at scales larger than about one- 
third the size of the shorter side. 

Providing that the image does not truncate the clus- 
ter emission in any direction, this shape problem can be 
treated as a classical edge effect. A simple solution is then 
to use the asymptotic background values to obtain a sim- 
ulated larger square image. From a practical point of view 
only a small truncation of the data occurs in our images 
along the shorter axis. Its impact on an azimuthal pro- 
file can therefore be considered as negligible and we thus 
decided to apply the above procedure to modify the ge- 
ometry of the EUVE images. We now give details on the 
required background estimation. 

Fig. |l| shows the rebinned image of a long observation 
of the Virgo cluster. The cluster is clearly visible in the 
central area of the field of view. Three rectangles are su- 



r\ and T2) 1 one towards each edge, as far away from the 
central source as possible. Both r\ and are assumed to 
contain only noise and to be representative of the noise 
level far from any celestial source. This is ensured because 
the centers of both rectangles lie at ~ 40 arcmin of the 
central source, a distance notably greater than the de- 
tection radius of the extended EUVE cluster emission, as 
seen with classical methods. The second order moment of 
the count distribution is used to estimate the level of the 
noise in each of these rectangles. The noise level in the 
large box is then computed as the mean of the noise levels 
in both rectangles. Once this has been done, each pixel 
outside the larger box is replaced by a random number 
drawn from a Poisson distribution, the variance of which 
is obtained by adding in quadrature the variances in each 
of the rectangles. The accuracy of the noise estimation is 
readily observable in the product image: if the noise is un- 
derestimated, the central box will obviously rise against 
the rest of the image, with its boundaries showing a clear 
discontinuity. On the other hand, if the noise has been 
overestimated, the outside parts of the new image will 
rise against the central box, again showing a discontinuity 
on the boundaries. 

When more than one EUVE observation of a cluster 
was available, we used the co-added images in order to in- 
crease the signal to noise ratio. Fig. || illustrates the case 
of Abell 1795. This image is typical of the co-added images 
we dealt with in the case of other clusters. The difference 
in roll-angle of the satellite in both observations is appar- 
ent. The central intersection of both observations obvi- 
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ously has the highest signal to noise ratio. It corresponds 
to a region indicated by the circle of radius 16 arcmin in 
Fig. ||, which we will take as the radial limit for our analy- 
sis. For each of the separate observations, we drew again a 
box closely following the outer detector shape. Both boxes 
are labelled (A,B,C,D) and (A',B',C',D') in Fig. §. In each 
of the boxes, we define two small rectangles, where the 
noise is evaluated. The noise level in a box is computed as 
the mean of the noise levels in both rectangles. Once this 
has been done for each box, each pixel outside the intersec- 
tion of the boxes (the central diamond) is again replaced 
by a random number drawn from a Poisson distribution, 
as discussed before. Once again, the under or overestima- 
tion of the noise will be clearly visible. An example of such 
an image is shown in Fig. 

2.3. The ROSAT data 

The ROSAT X-ray data for each cluster were obtained 
from the HEASARC web archive]]. Owing to the higher 
sensitivity, we preferred PSPC to HRI data and pointed 
observations to all-sky survey scanning mode observa- 
tions. When multiple PSPC observations were available, 
we always used the one with the smallest offset radius 
compared to the center of the cluster, and then the longest 
one. 

S. Snowden's Extended Source Analysis Software 
(ESASQ) was used to perform the data reduction, i.e., re- 
jection of high-background times, modelling of the differ- 
ent background types, energy-dependent background sub- 
traction, exposure and vignetting corrections (see Snow- 
den et al. 1994). We adopted a conservative value of 
170 cts s _1 for the maximum Average Master Veto rate 
allowed. We then carefully examined the light curves of 
the total counts in the entire image per energy band, and 
checked that no short time scale glitches were present. 
We limited our analysis to the [0.5,2.0] keV energy band 
(bands R4-R7 as defined by Snowden). This prelimi- 
nary reduction produced a surface brightness image with 
512 x 512 pixels, 15 arcsec per pixel (roughly the FWHM 
of the PSPC PSF at 1 keV), as well as an exposure im- 
age for each of the bands and a background image. The 
image was then divided by the mean exposure map. We 
thus obtained a 512 x 512 pixel image, in units of counts 
s^ 1 pixel -2 , which was then rebinned to 18. 5x 18.5 arcsec 2 
pixels to match the EUV superpixel size. 

The ROSAT images do not present the same problems 
as the EUVE ones, because in this case, the geometry 
is circular around the central cluster (all our observations 
were pointed). The only geometrical constraint arises from 
the presence at an offset angle of approximately 20 arcmin 
of the structure support of the telescope. However, since 
this value is smaller than the EUVE limiting radius, it 

1 http://heasarc.gsfc.nasa.gov/ 

2 ftp:/ /legacy. gsfc.nasa.gov/rosat/software/fortran/sxrb 



does not place any further constraint here. We therefore 
just extracted a square image inscribed in the 20 arcmin 
radius circle. 

2.4- Wavelet analysis and reconstruction 

The comparison between the EUVE observed surface 
brightness profiles of clusters and their ROSAT X-ray pre- 
dicted counterparts is the key step in our analysis. How- 
ever a direct comparison of the X-ray and EUVE data 
would provide misleading results due to the noise in the 
raw data. To remove most of this noise while preserving 
information on the astronomical emission at all spatial 
scales of interest we decided to apply part of the wavelet 
vision model described in detail in Rue & Bijaoui (1997). 

Let us hereafter briefly summarize the main steps of 
the image processing involved. First, a discrete wavelet 
transform of the image Iq is performed using the so-called 
a trous implementation, which is nothing more than the 
classical multiresolution algorithm of Mallat (1989) but 
without any decimations since "holes" are introduced in 
the required convolutions. The initial sampling is kept 
for each scale; this insures shift-invariance and allows 
straightforward interscale computations. The scaling func- 
tion leading through a dilation factor a to the involved sets 
of low-pass and band-pass (wavelets) filters is a cubic B- 
spline (defined in the Fourier space as sinc 4 (w)), very well 
suited for enlightening gaussian-like patterns of various 
sizes in an image. Owing to the properties of the under- 
lying convolutions a stationary white noise signal leads to 
a statistical distribution of the wavelet coefficients with a 
zero mean value at each computed scale aj and a stan- 
dard deviation decreasing as the spatial scale increases. 
Structures in the transform which are related to objects 
in direct space can then be identified and located by se- 
lecting high enough wavelet coefficients at each scale. 

Removing noise from the data without any loss in res- 
olution can thus be achieved quite easily in wavelet space 
by keeping only those coefficient values with a probability 
to be due to a chance fluctuation of the noise process lower 
than a chosen decision level. This second step of the noise 
removal process requires the knowledge of the probability 
density function of the coefficients at scale a\ = 1 pixel 
(here 18.5 arcsec) for the noise process and of the de- 
creasing law of the rms scatter a, the asymptotic value of 
which appears to be 0.5 for a bidimensional image Iq and 
a dyadic scale analysis (aj = 2 J_1 , j e Z) : <jj = a\jV^ x . 
In case of Gaussian noise, the probability density function 
is easily obtained from a histogram of the wavelet coeffi- 
cients providing that this first scale of one pixel in size is 
noise dominated. Due to the stationarity of the Gaussian 
process, the selection of the so-called statistically signif- 
icant coefficients Vj at scale aj can then be done using 
the same k-aj threshold. We end up in this way with a 
subset of coefficients at different scales, the pixels of 
which define what is called hereafter the multiscale signif- 
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icant mask Mq. In case of Poisson noise stationarity can 
be obtained by applying an Anscombe (1948) transform 
on the data : f(x) = 2\Jx + 0.375. When the mean pa- 
rameter fi of the Poisson process is greater than about 5, 
this transform indeed leads to a new process with a very 
small bias in the expectation value ( < ^ X j^ f 7)^^ < 
and a nearly perfect variance stabilization. The same de- 
cision scheme as in the Gaussian case can therefore be 
used for defining the multiscale significant mask from the 
wavelet transform of this image with a constant variance. 
Only second order errors in the wavelet coefficients would 
arise at the very first scales in areas with very low counts. 
This is not a matter of concern for the present application 
since the computation of a correct azimuthal brightness 
profile necessitates only a proper account of the medium 
and large-scale features in these faint regions. 

The third and final step is to restore a positive image I 
where the noise is strongly reduced from the subset v 1 ^ 1 of 
the significant wavelet coefficients of 1^ = Iq- This is an 
inverse problem which has in general an infinity of solu- 
tions and the choice among these solutions must be done 
with respect to a regularization constraint. To do so, an 
iterative conjugate gradient algorithm is used. This algo- 
rithm minimizes the distance ||z/ )— where z/ n ) is the 
restriction to Mq of the wavelet transform of the solution 
at step n. Stating A the operator associated to the 
wavelet transform followed by the projection onto Mo and 
A the joint operator, = A(z/™)) is obtained from 

the previous solution by : = /("■) + a^Ir where 

the residual imag e 4 n) is defined as A(^ n) ) + 
with Vr = v^ — A(/(™)) ; a and (3 are convergence pa- 
rameters. Convergence ( ^"jjlj^i] ^ < ^ s reacne d after 
very few steps, leading to an image / without any discon- 
tinuities despite the hard thresholding performed in the 
wavelet space. Note also that any flat instrumental or as- 
tronomical component or gradient in the initial image is 
automatically removed in image / taking benefit of the 
mathematical properties of the wavelet function related 
to our scaling function. 

We applied such an image processing to both our 
EUVE and ROSAT square images considering seven or 
eight spatial scales for the wavelet transform according to 
the size of the image. Pure Poisson noise was considered 
in both cases. This is not strictly the case for the ROSAT 
data because of the preprocessing of the raw counts (ex- 
posure correction, modelled background subtraction) , but 
this approximation appears to be good enough for detect- 
ing correctly all the large-scale components we are inter- 
ested in. Significant wavelet coefficients were thus selected 
according to a probability of false detection of 4cr, 3.5a and 
3cr for scales 1 (18.5 arcsec), 2 (37 arcsec) and 3 (74 arcsec) 
and above respectively. The rather high levels for the first 
two scales automatically eliminate most of the small-scale 
structures if any (these might distort the overall shape of 
the brightness profile or might have arisen from an error 



in the probability density function of the noise or from a 
loss in control over the bias in the Anscombe transform in 
regions with very few counts). We have checked in several 
cases that the profiles do not change significantly when 
the thresholding for scales 1, 2 and 3 is modified. 

2.5. Radial profiles 

After wavelet reconstructing both EUVE and ROSAT 
images for each cluster, we proceed to compare the de- 
crease of the surface brightness in the two observational 
bands. Profiles were thus derived from the 2D recon- 
structed images within concentric elliptical rings using 
the STSDAS. ANALYSIS. ISOPHOTE.ELLIPSE task in 
IRAF. The cluster ellipticities and major axis position 
angles were estimated from the X-ray images, which have 
better signal to noise ratio, except for Coma because it 
fills the 20 arcmin PSPC ring and its cllipticity had to be 
derived from the EUVE reconstructed image. The ellip- 
ticities and major axis position angles were then fixed to 
be the same in the EUV and X-rays, in order to measure 
intensities in the same spatial regions. 

In order to see how the instrument point spread func- 
tions (PSF) could influence our results, we also analyzed 
a point source in a similar way: the star 7Tau, for which 
both EUV and X-ray data were available. The EUV and 
X-ray exposure times for this star were 89924s and 8547s 
respectively, and the ROSAT PSPC image of 7Tau was di- 
vided by its corresponding exposure map, as for clusters. 
The EUV and X-ray images were rebinned to a pizel size of 
18.5x18.5 arcsec 2 as previously, and the EUV and X-ray 
profiles were then derived in a similar way. Note that this 
point source is about 15 arcmin offset from the detector 
center in the ROSAT image, at a radius where the PSF 
of PSPC is already degraded. Therefore, the PSF effect in 
X-rays will in fact always be smaller or equal than what 
we measure for 7 Tau. This is not the case for EUVE, 
where the point source is close to the detector center. 

The problem of a possible non-flatness of the back- 
ground was already addressed in the case of EUVE by 
Lieu et al. (1999c). A comparable wavelet analysis was per- 
formed on a blank field totalling 85 ksec exposure time, 
and the backgroud was found to vary by at most 10% 
within the central circle of 16 arcmin radius defined above 
taken into account here (see Fig. 1 by Lieu et al. 1999c). 
The mean value of another blank field observed for 21 ksec 
is found to be 1.21 counts/s/pixel 2 , excluding the outer- 
most regions (that is in a rectangle similar to ABCD from 
Fig. [l]), which is comparable to the mean background value 
in Abell 1795 (1.19 counts/s/pixel 2 ). 

In the case of the X-ray images, they were divided by 
their corresponding exposure maps, and therefore the non- 
flatness of the ROSAT PSPC detector should be corrected 
for. We are therefore confident that the non-flatness of the 
detector backgrounds cannot change influence our results 
by more than a few %. 
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Error bars on all profiles are those directly given by the 
ELLIPSE task. However, since we were concerned that 
they might be underestimated, we also computed them 
from the Poisson fluctuation in the number of counts in 
each ring ; these error bars were found to be commensu- 
rate with those given by ELLIPSE (the rms fluctuation of 
the pixel intensity within each concentric ellipse). For re- 
constructed images, they are typically smaller than ±0.01 
in logarithmic scale, i.e. too small to be clearly visible on 
Figs. | to |). 



3. Results 

3.1. EUV and X-ray profiles 

For co-added images, we do not expect to detect any struc- 
ture outside a certain radius since the outer part of the 
image is purely noise. The central circle in Fig. |^ shows 
a conservative estimate of this radius, i.e. the half width 
of the EUVE detector ~16 arcmin, which we will take as 
the limit of our analysis. 



Radius (kpc.) 
500 1000 1500 




5 10 15 

Radius (arcmin) 



Fig. 4. EUV and X-ray profiles of Abell 1795, obtained 
directly from the raw data (points) and after wavelet anal- 
ysis and reconstruction (points with connecting lines). 
As in the four following figures, the EUVE data are 
drawn with triangles (bottom curves) and the X-ray data 
with filled squares (top curves); the radius is expressed 
in arcminutes (bottom) and kpc (top); intensities are in 
counts/s/pixel 2 after background subtraction (1 pixel is 
18.5x18.5 arcsec 2 ). The error bars are drawn, but for 
the profiles derived from the wavelet reconstructed images 
they are too small to be clearly visible (typically smaller 
than ±0.01 in logarithmic scale for Figs. | to |). The EUV 
(crosses) and ROSAT PSPC (empty squares) PSFs drawn 
from the images of the point source 7Tau, shifted by an 
arbitrary value to make the figure clearer, are shown at 
the bottom for comparison. 



We first derived the profiles of the EUVE and X-ray 
images of all clusters before and after wavelet analysis and 
reconstruction, as a test to check that our method con- 
serves the overall profile shape and intensity level. The 
corresponding profiles are shown in Fig. [|for Abell 1795. 
Note that all the profiles drawn in Figs. |]to ||are obtained 
after subtracting the background. For the raw data, this 
background is estimated as a mean value far from the 
edges of the detector; for wavelet reconstructed images, 
the background is modelled during the wavelet analysis 
itself and subtracted. There is good consistency between 
the two EUVE curves on the one side and between the two 
X-ray curves on the other side, implying that the wavelet 
analysis and reconstruction applied here does not affect 
fluxes, but simply provides a means of detecting cluster 
emission out to a somewhat larger radius (as in the case 
of the EUVE data). The small feature apparent in the X- 
ray profile at ^6 arcmin is due to a Seyfert galaxy super- 
imposed on the cluster; it was removed from the wavelet 
reconstructed image and therefore does not appear on the 
corresponding profile. Note however that Abell 1795 is the 
only cluster for which the raw and wavelet reconstructed 
X-ray profiles do not coincide exactly (in the other clus- 
ters they can hardly be distinguished). Although we have 
tried various thresholding values for the small scales and 
various background subtractions of the raw data to see 
if they could coincide better, we have not been able to 
superimpose them exactly. 

The EUV and ROSAT PSPC PSFs drawn from the 
images of a point source (see Sect. 2.5.) for comparison are 
also plotted in Fig. || They indicate that the PSFs have 
no influence on the profile shapes for radii larger than 
about larcmin and 2.5 arcmin in the EUV and X-rays 
respectively. Therefore any features found in the profiles 
within regions smaller than these radii are probably not 
real. 

The EUV and X-ray radii of Abell 1795 are at least 
16 arcmin, corresponding to a physical extent of about 
1700 kpc. This extent is notably larger than the previously 
found value of 10 arcmin (Mittaz et al. 1998, Bonamente 
et al. 2001b). 

The EUV and X-ray radii of Abell 2199 extend to at 
least 16 arcmin, or 820 kpc (Fig. ||), comparable to the 
value of 15 arcmin measured by Kaastra et al. (1999) both 
with EUVE and Beppo-Sax. 

In Abell 4059 (Fig. |), the EUV and X-ray emissions 
also reach 16 arcmin or 1300 kpc, much further than the 4- 
5 arcmin reported by Berghofer et al. (2000b). The slopes 
of the profiles in the EUV and X-rays differ notably, at 
least in the inner 8 or 10 arcmin. The difference in the 
central regions between the EUV profiles derived from the 
raw and wavelet reconstructed images may be due to the 
fact that structures only present at small scales were elim- 
inated, and therefore signal may be missing in the very 
center. This cannot be due to the shape of the EUVE 
PSF, since it is flat for radii larger than about 1 arcmin 
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Fig. 5. Same as Fig. | for Abell 2199. 



Radius (kpc) 
500 1000 



Abell 4059 



■■•'I... 



/:::;7''":"'-«-'-:-,........ : ..:.; :: ; 



o -6 - 

5 



i i i i I i i i , I 

5 10 

Radius (arcmin) 



Fig. 7. Same as Fig. [| for Coma. 
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Fig. 6. Same as Fig. | for Abell 4059. 



Fig. 8. Same as Fig. [|for Virgo. 



while this discrepancy is observed up to a radius of about 
3 arcmin. Such a behaviour is not observed in Abell 1795, 
Abell 2199 or Coma, probably because they are brighter. 

The EUV and X-ray radii in Coma reach at least 16 
arcmin (Fig. 0), corresponding to 630 kpc and slightly 
larger than reported by Bowyer et al. (1999). The two 
profiles appear quite similar at least up to a radius of 
about 10 arcmin. 

Virgo is the most nearby cluster in our sample. EUV 
and X-ray emissions reach at least 16 arcmin (Fig. ||), that 
is slightly larger than detected by Lieu et al. (1996a) and 
Bcrghofcr et al. (2000a), but smaller than the 20 arcmin 
claimed by Bonamente et al. (2001b); however, this extent 
corresponds to the very small physical value of ~ 80 kpc, 
due to the proximity of Virgo. Both profiles appear quite 
parallel except in the inner 4 arcmin; note however that 
since the X-ray PSF is not flat for radii smaller than 2.5 
arcmin, the difference between the EUV and X-ray profiles 
in the inner 4 arcmin may not be real. 



In order to quantify the probability for the EUV and 
X-ray profiles to be similar, we performed a Kolmogorov- 
Smirnov (hereafter K-S) test on the EUV and X-ray pro- 
files drawn from the wavelet reconstructed images, both in 
concentric elliptical annuli, and in ellipses with increasing 
radius, after normalizing the EUV and X-ray profiles to 
the same innermost pixel value. 

In all clusters, the probability that the two profiles are 
issued from the same parent population varies with radius: 
it is a few % in the very center (1.5 arcmin radius), and 
less than 2% in concentric annuli of width 1.5 arcmin (i.e. 
for radii between 1.5 and 3.0, 3.0 and 4.5 arcmin etc.). 
If we consider ellipses of radii 3.0, 4.5, 6.0 arcmin etc. 
instead, the probability that the two profiles are issued 
from the same parent population is smaller than 0.1% in 
all cases. Therefore even if the EUV and X-ray intensity 
profiles may appear similar to the eye, as for example in 
Coma, it is clear that statistically they are not. 
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3.2. EUV to X-ray intensity ratios 

In order to compare the EUV to X-ray profiles, we com- 
puted the EUV to X-ray intensity ratios in the same con- 
centric ellipses; these are shown in Figs. || to [l3| The ratio 
of the EUV to X-ray PSFs was also derived (see Fig. ^). 
As expected from the individual EUV and X-ray PSFs, 
this ratio is only constant for radii larger than 2.5 arcmin, 
while below it shows a strong dip. Therefore the EUV to 
X-ray ratios below this radius are not reliable. Such ratios 
have to be compared to those expected in the hypothesis 
that the EUV emission is the low energy tail of the hot 
ICM emission. We therefore employed spatially resolved 
hot ICM models for each cluster, in order to account for 
the effects of temperature and abundance gradients on 
the predicted ratios. Spectra of all cluster regions were 
modelled with a thin-plasma emission code (MEKAL in 
XSPEC, Mewe et al. 1985, 1986, Kaastra 1992) modified 
by photoelectric absorption (WABS in XSPEC, Morrison 
& McCammon 1983). In detail, the hot ICM in Abcll 1795 
was modelled with kT increasing from 3.5 keV at the cen- 
ter to 6 keV at radii > 3 arcmin, and abundances in the 
0.3-0.5 solar range (Briel & Henry 1996); Abell 2199 was 
modelled with kT in the 2.5-4 keV range, and abundances 
decreasing from ^0.7 solar at the center to 0.3 solar for 
radii > 3 arcmin (Siddiqui et al. 1998); Abell 4059 was 
modelled with kT increasing from 2 keV at the center to 
5 keV at radii > 3 arcmin, and abundances in the 0.3-0.6 
solar range (Hwang & Sarazin 1998); Coma was considered 
as isothermal with kT=8.2 keV and abundances 0.21 solar 
(Hughes et al. 1993); finally, the hot ICM parameters of 
Virgo were directly fit to our PSPC data, with kT in the 
1.2-3 keV range, and metal abundances of 0.4-0.6 solar. 
Using these models, we were able to derive the expected 
count rates in the EUV DS band and in the PSPC 1/4 
keV band, taking into account the effective area of the in- 
struments and the performances at all detector positions. 
The points corresponding to these calculated count ratios 
are shown in the following figures. 

In Abell 1795, the EUV to X-ray intensity ratio (Fig. §) 
is smaller than 1 in the center (for radii smaller than 1.5 
arcmin), but this is most probably due to the shape of the 
ratio of the EUV to X-ray PSFs. It then increases with 
radius to reach a roughly constant value of 4. The shape 
of this curve is consistent with that given by Bonamente 
et al. (2001b), while the value of 4 that we find in the 
outer regions is intermediate between that of 6 given by 
Mittaz et al. (1998) and that of 2 given by Bonamente et 
al. (2001b). Error bars are typically smaller than ±0.01 
in logarithmic scale, and are therefore difficult to see on 
the graphs. We therefore confirm the strong soft excess 
found in this cluster by these two groups of authors, and 
disagree with Bowyer et al. (1999) who claimed that no 
excess EUV emission was seen in this cluster. Note that 
the small disagreement between the raw and wavelet re- 
constructed X-ray image profiles (Fig. flj) cannot be at the 
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Fig. 9. Ratio of the EUV to X-ray emission in Abell 1795 
in logarithmic scale. As in all following figures, squares 
show the measured EUV to X-ray flux ratios (error bars 
are difficult to see on the graphs since they are typically 
smaller than ±0.01 in logarithmic scale, as for the other 
clusters in Figs. |ili|-|l3|). Triangles indicate the EUV to X- 
ray flux ratios in various concentric rings, expected from 
the hypothesis that the EUV emission is the tail of the 
bremsstrahlung emission accounting for the observed X- 
rays (see text). The ratio of the EUV to X-ray PSFs is also 
drawn as empty squares, shifted by an arbitrary value to 
make the figure clearer; to guide the eye, the mean value 
of this ratio is indicated with a dashed line. 
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Fig. 10. Same as Fig. | for Abell 2199. 



origin of the EUV excess that we observe, since on the con- 
trary the raw X-ray flux is higher that the reconstructed 
one, and therefore using the raw X-ray data would lead to 
an even larger EUV excess. 

The EUV to X-ray intensity ratio in Abell 2199 is 
smaller than 1 for radii below 5 arcmin, but as for Abell 
1795 this is most probably due to the shape of the ratio of 



Florence Durret et al.: Excess EUV emission in five clusters 



9 



Radius (kpc) 
500 



Abell 4059 



5 10 
Radius (arcmin^ 



Radius (kpc) 

20 40 60 80 

i i | i i i | i i i | i i i r 

Virgo 



5 10 
Radius farcmin) 



Fig. 11. Same as Fig. | for Abell 4059. 
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Fig. 12. Same as Fig. |^ for Coma. 

the EUV to X-ray PSFs. It then increases and reaches a 
value of about 3.7 between radii of 10 and 16 arcmin. We 
therefore agree with the soft excess observed in this clus- 
ter from various sets of data by Kaastra et al. (1999), and 
here also disagree with Bowyer et al. (1999) who found no 
excess EUV emission in this cluster. 

The EUV to X-ray intensity ratio in Abell 4059 is 
smaller than 1 for radii below 2 arcmin (Fig. [ll]), in agree- 
ment with Berghofer et al. (2000b), but again this is most 
probably due to the shape of the ratio of the EUV to 
X-ray PSFs; for increasing radii, it strongly increases and 
reaches a value of almost 20 between radii of 12 and 16 ar- 
cmin, in total disagreement with Berghofer et al. (2000b) , 
who find no EUV emission beyond 4 arcmin and no EUV 
excess beyond 2 arcmin. 

Contrary to what happens in the three previously stud- 
ied clusters, the EUV to X-ray intensity ratio in Coma 
appears rather flat up to a radius of at least 10 arcmin 
(Fig. |lj), but as shown above the two distributions are 



Fig. 13. Same as Fig. ^ for Virgo. 

not similar. The fact that no dip is seen within the central 
few arcminutes, as would be expected from the shape of 
the ratio of the EUV to X-ray PSFs could indicate that 
the EUV excess in the center is even stronger than seen 
in Fig. [l2|. There is an EUV excess by a factor of ~ 1.4 
over a thermal bresmsstrahlung tail. 

The EUV to X-ray emission ratio in Virgo is in 
agreement with the EUV excesses previously reported by 
Berghofer et al. (2000a) and Bonamcntc et al. (2001b), 
therefore implying that here also the EUV emission is not 
a thermal bremsstrahlung tail. The rise and fall of the 
EUV to X-ray ratio for radii increasing between 0.3 and 
2.5 arcmin (Fig. |l3| ) is rather surprising since it varies in 
an opposite way to the ratio of the EUV to X-ray PSFs; 
however images with higher spatial resolution are obvi- 
ously necessary to push the analysis further at these radii. 

The superimposed EUV to X-ray intensity ratios for 
all the clusters of our sample are plotted in Figs. |lj and[l5| 
for radii expressed in arcminutes and in kpc respectively. 
Although the shapes of these curves in the central 2.5 
arcmin most probably reflect that of the ratio of the EUV 
to X-ray PSFs (Fig. |lj), we can see that at least beyond 
this radius all the clusters do not have the same behaviour. 
This is particularly true when radii are expressed in kpc, 
and certainly must be linked to a difference in the emission 
mechanisms of the EUV emission in these clusters. Note 
that the radial extent of the EUV and X-ray emission in 
Virgo is comparable to the other clusters when the radius 
is expressed in arcminutes, but becomes extremely small 
in physical distance units, due to the very small redshift 
of the cluster: slightly more than 80 kpc. 

4. Discussion and conclusions 

In this paper, we have shown unambigously the existence 
of a EUV excess in all five clusters of our sample. 

In the first three clusters (Abell 1795, Abell 2199 and 
Abell 4059) , the EUV to X-ray intensity ratios have com- 
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Fig. 14. EUV and X-ray profiles for the five clusters in our 
sample, with the radius expressed in arcminutes. The sym- 
bols are the following: Abell 1795: filled triangles, Abell 
2199: empty squares, Abell 4059: empty circles, Coma: 
empty triangles, Virgo: filled squares. Error bars are omit- 
ted for clarity. The dashed line shows the ratio of the EUV 
to X-ray PSFs. 
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Fig. 15. EUV and X-ray profiles for the five clusters in 
our sample, with the radius expressed in kpc. 



parable behaviours: they show a possible deficiency of 
EUV emission over a bremsstrahlung tail in the very cen- 
tral regions (but this may be an effect of the ratio of the 
EUV to X-ray PSFs), and a EUV excess beyond a few ar- 
cmin. As suggested by several authors (Bowyer et al. 1999; 
Lieu et al. 2000), the first of these features, if real, can be 
interpreted as due to excess absorption within the cluster 
core due to the fact that in the cooler central regions some 
metals are not fully ionized and these ions absorb part of 
the soft X-ray flux. The presence of such excess absorption 



would agree with the previously claimed existence of cool- 
ing flows in these three clusters. However, XMM-Ncwton 
has detected much weaker emission lines than expected 
from bona fide isobaric cooling flow models, implying that 
there is significantly less "cool" gas than predicted (e.g. 
Kaastra et al. 2001). Therefore it is no longer straightfor- 
ward to interpret excess absorption at the center of these 
clusters as due to the presence of cooler gas. 

One may however note that the only cluster (Coma) 
for which no EUV dip is seen in the very center is both 
the hottest one by far and the only one with no "cooling 
flow" whatsoever. Since the ratio of the EUV to X-ray 
PSFs shows a dip for radii smaller than about 2.5 arcmin, 
the absence of such a dip in Coma suggests that in fact 
there is a significant EUV excess in the central regions of 
Coma. 

Virgo shows a different behaviour, since the EUV ex- 
cess is the strongest in the very center, within radii smaller 
than about 4 arcmin (corresponding to a physical size of 
~20-30 kpc, that is roughly the size of M87), then remains 
roughly constant with radius between 4 and 7 arcmin, and 
finally increases with a shallow slope for radii larger than 
7-8 arcmin. As for Coma, the presence of a bump in the 
central regions suggests that the EUV excess is in fact 
very stronge in this zone. For these last two clusters, im- 
ages with higher spatial resolution are obviously necessary 
to analyze the EUV excess at these small radii. 

The mere presence of an EUV excess in all five clus- 
ters of our sample indicates that a mechanism other than 
bremsstrahlung is needed to account for this emission. Our 
data does not allow us to discriminate between the various 
mechanisms proposed in the literature for the soft excess. 
However, in view of the most recent results obtained with 
XMM-Newton suggesting that there is much less warm 
gas in the central regions of clusters than previously be- 
lieved, it seems likely that this EUV excess is probably of 
non thermal origin. 
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